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High frequency QPOs: nonlinear oscillations in strong gravity 
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Abstract. Some of the more promising ideas about the origin of the high frequency variability (kHz QPOs) in the observed 
X-ray emissions of low-mass X-ray binaries are contrasted with less promising ones. 
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1. Setting the stage 

1.1. Remembering the key crossroads 

If I am personally convinced that the high(est)-frequency 
QPOs in black hole and neutron-star sources have a common 
origin, it is not only because this assumption has served us 
well in the past five years, but also because I have taken the 
other road, and I have seen that it leads nowhere. I would like 
to begin by describing my journey, so that others may decide 
whether or not the direction we now take is reasonable. 

1.2. The accretion-gap paradigm for neutron stars 

The discovery twenty years ago of the first millisecond (ra- 
dio) pulsar immediately brought home the idea that low- 
mass X-ray binaries (LMXBs) may contain rapidly rotating 
neutron stars endowed with a magnetic field far below the 
theretofore canonical value of 10^^ G. These were hard-to- 
imagine times when the timing capability of the the best X- 
ray instrument (GINGA) was below ~ 200 Hz, and (EX- 
OSAT) spectra were only reported up to 30 keV. 

In our first paper on LMXBs, formally devoted to in- 
vestigating spin-up by accretion, we examined the general- 
relativistic orbits around neutron stars with dynamically 
unimportant magnetic fields. We have pointed out that for 
many viable equations of state (of matter at supranuclear 
densities) the accretion disk could terminate close to the 
marginally stable orbit (ISCO) well outside the stellar sur- 
face (Kluzniak and Wagoner 1985) and dubbed this situation 
the "accretion-gap regime." In other words, the accretion disk 
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of a neutron star in a typical LMXB could well be very sim- 
ilar to the accretion disk around a black hole, with the inner 
edge determined by GR effects of strong gravity (with frame 
dragging and all). 

The next paper in the short series, prepared as part of my 
thesis work (Kluzniak 1987), demonstrated that the spectrum 
of the emission expected in the accretion-gap regime, when 
the freely falling matter which left the inner edge of the disk 
hits the equatorial accretion belt, is a rather hard power-law 
extending up to ~ 10^ keV (Kluzniak and Wilson 1991). The 
discovery several years later of such emissions from some 
neutron star sources can be taken as possibly confirming the 
presence of the ISCO. However, non-thermal spectra do not 
uniquely point to the geometry or process of emission, and so 
we also searched for timing signatures of the accretion-gap 
regime or, as it is put today, of a key prediction of Einstein's 
theory of gravitation (GR). 

We were looking not for a tiny quantitative departure 
from Newtonian gravity, but for a major difference. The time 
of flight of matter through the gap would be one such signa- 
ture of a qualitative departure from Newtonian accretion ge- 
ometry. Accordingly, we computed the trajectory and the fate 
of a clump of matter orbiting in the accretion disk, crossing 
the sonic point within the inner edge, and hitting the stellar 
surface (Kluzniak, Michelson and Wagoner 1990). We pro- 
posed this as a method of measuring the mass of the neutron 
star, or (if the mass were known) of testing GR. In particu- 
lar, interpreting the 0.2 Hz QPO observed in a certain X-ray 
pulsar (Angelini, Stella and Parmar 1989) as (an unknown 
process) occurring at the orbital frequency close to the inner 
boundary of the accretion disk, which is ^ 10'^ stellar radii 
above the surface of this strongly magnetized accreting neu- 
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tron star, we pointed out that in a LMXB neutron star in the 
accretion gap regime the same QPO should occur at a fre- 
quency which a) is in the kHz range, b) inversely depends on 
the mass M of the neutron star, and c) depends on the angular 
momentum j of the neutron star — the frequency at the ISCO 
being 2.2kHz {Mq/M){1 + 0.749j). It is in the context of 
the accretion-gap paradigm that our efforts to understand kHz 
QPOs should be understood. 

1.3. Two kHz QPOs in neutron stars, one in black holes? 

The discovery in neutron stars of two kHz quasi-periodicities 
in observations of neutron stars by the first instrument 
(RXTE) that had the capabiUty of observing such rapid vari- 
ability did not seem very surprising in light of the Kluzniak, 
Michelson and Wagoner 1990 prediction discussed above. 
Neither did the lack of a coherent pulsation associated with 
the spin of the neutron star — after all, to obtain high orbital 
frequencies, we had to assume that the stellar magnetic field 
was too weak to channel the flow. The frequencies came out 
just right, if the neutron stars were a bit on the heavy side 
of the allowed mass range in an LMXB transferring mass for 
~ 10* y. 

Soon a high frequency QPO was reported also from a 
black hole source. This all seemed to make perfect sense. 
Black hole acretion disks were just like neutron star acre- 
tion disks, the only difference being in the size of their in- 
ner gap (proportional to mass). Hence, one QPO at the inner 
edge of the disk was expected, probably associated with an 
orbiting clump. On top of it, matter falling onto the equato- 
rial accretion belt differentially slipping on the neutron star 
surface — i.e., rotating faster than the neutron star, albeit at 
variable speeds — could support a structure (perhaps a vortex 
like the great red spot on Jupiter, although the referee pre- 
vented this specific suggestion from appearing in print) giv- 
ing rise to the second quasi-period. No such second QPO was 
expected from a black hole, where matter fell down the hori- 
zon. At any rate, identifying the observed frequency with that 
in the ISCO one could determine the mass of the neutron star 
or black hole (Kluzniak 1998), a favorite pastime for many 
authors in the subsequent years. 

A model very similar to this (minus the strong gravity 
part), had been proposed much earUer by Paczynski to ex- 
plain the white dwarf QPOs, and their variable frequencies, 
so it all made perfect sense in the kHz QPO context, where 
both frequencies were known to increase with the luminosity 
of the source — the idea in the model being that at higher ac- 
cretion rates the sonic point moves in closer to the stellar sur- 
face, and that the top of the equatorial accretion belt moves 
faster as the accretion rate of angular momentum increases 
(see Kluzniak, Michelson and Wagoner 1990 for references 
on the position of the sonic point). 

However, when I tried to make the model of the frequency 
variation in kHz QPOs agree quantitatively with the data, I 
failed. And I could make no sense of any other published 
model based on the assumption that the presence of the sec- 
ond kHz QPO is a direct signature of the (spinning) stellar 
surface. A timely warning fromDimitrios Psaltis that the final 



word on the number of QPOs in black holes has not yet been 
said, as they are very weak at high frequencies, has helped 
me to finally abandon the idea that only one of the high fre- 
quency QPOs arises in the accretion disk. If one arises in the 
accretion disk, why not both? 

1.4. Wrong turn 

Much initial work on netron-star QPOs was informed by ear- 
Uer models of similar phenomena in white dwarf systems 
(CVs), but the assumption that the kHz QPOs were caused 
by orbiting clumps (inhomogeneities in the accretion flow) 
turned out to be wrong. It had been known that clumps in 
the differentially rotating accretion disk would be sheared out 
into nearly axisymmetric structures within a few orbits (Bath, 
Evans and Papaloizou 1974) — i.e., a clump is short lived — 
and it has simply been overlooked that the kHz QPOs were 
reported to be very narrow features already in the discovery 
papers. Today we know that the high quality factor, Q ~ 200 
of the lower kHz QPOs in neutron-stars excludes all clump- 
based models, as well as those involving longer-lived vortices 
drifting radially inwards with the mean accretion flow (Barret 
et al. 2005). Thus, the question of relative merits of the "beat 
frequency model" and the "relativistic precession model" is a 
moot point. 

1.5. Linear diskoseismology 

All the while, fundamental work was being done on the nor- 
mal modes of accretion disks. It turns out, that in GR g-modes 
are trapped near the maximum of the radial epicychc fre- 
quency (Okazaki, Kato, Fukue 1987), an effect clearly absent 
in Newtonian 1 /r gravity that requires the epicyclic frequen- 
cies to be equal to the Keplerian orbital value ^jGMjr'^, 
itself a monotonic function of the radius. After nearly two 
decades of concerted effort (Wagoner 1999, Kato 2000) sev- 
eral modes have been identified, which were promising in the 
context of black hole HFPOs. Because this is such a difficult 
problem aU work had been carried out in the linear regime. 

One conclusion, more clearly seen a posteriori, is that no 
special ratio of frequencies for any of the preferred low-lying 
modes is predicted. More precisely, over and above the 1/M 
scaling in GR, the frequencies of the various modes are a 
function of the spin of the central compact object, and so is 
the ratio of frequencies, at least for modes which are not har- 
monics of the same fundamental. Thus, judiciously choos- 
ing the black hole spin one can acconmiodate nearly any fre- 
quency ratio of, e.g., the fundamental g-mode and c-mode. 
However, no two black holes (in a limited sample) are ex- 
pected to have the same spin, hence any such ratio must be 
accidental and limited to one or two objects at most. 

Another conclusion was that the frequencies have only a 
slight dependence on the expected variations in the disk struc- 
ture. Therefore, disk modes were not thought to be directly 
responsible for modifying the X-ray flux in neutron-star sys- 
tems, where the QPO frequencies vary quite a bit in any given 
source. 
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These two conclusions need no longer hold in non-linear 
diskoseismology. The equations of hydrodynanoics are intrin- 
sically non-Unear and there is urgent need to carry out an 
analysis of disk modes in the non-linear regime. We have rea- 
sons to beUeve that the properties of the high frequency QPOs 
in both black hole and neutron-star systems are a manifesta- 
tion of non-Hnear coupling between certain oscillations of the 
accretion disk. 

2. Two characteristic variable frequencies 

The two kHz QPOs in neutron stars have reproducible prop- 
erties in each source, including characteristic frequency val- 
ues, and yet the two frequencies are not fixed — ^they vary. Il 
is easy enough for the accreting fluid to display characteris- 
tic frequencies. Even in Newtonian gravity, where the orbital 
frequency is scale-free, there is a maximum frequency corre- 
sponding to that in the inner disk — simply because there is a 
characteristic radius, for example the steUar radius. In GR or- 
bital motion characteristic frequencies are built in, for exam- 
ple the maximum radial epicyclic frequency, and the frequen- 
cies of the trapped modes are a reflection of this property of 
strong gravity. However, all these characteristic frequencies 
are fixed for a given space-time metric (i.e., a given neutron 
star or black hole). 

It is the variation of frequencies that is highly suggestive 
of non-Unear dynamics. Frequencies that vary in a character- 
istic range are a well studied feature of non-linear oscillators. 
For these reason, we have suggested that the two kHz QPOs 
in neutron stars are a manifestation of a nonUnear resonance 
in an accretion disk placed in strong gravity (Kluzniak and 
Abramowicz 2001a,b). Two corrolaries immediately follow, 
and within a few months both were directly verified by ob- 
servers unaware of our hypothesis. The first is that there are 
preferred rational ratios of frequencies, such as 1:2 or 2:3. 
The second is that high frequency QPOs in black holes also 
come in pairs. 

This triumph of simple physics invites more work to be 
done. In paricular, it would be good to identify the two non- 
linear oscillators, and it is necessary to convincingly identify 
the similarities and the differences of black hole and neutron- 
star QPOs. This is the subject of the current workshop. 

3. The smoking gun of non-Unear oscillations 

The twin high-frequency QPOs need not have identically the 
same propeties in black holes and neutron- stars. The differ- 
ences, if any, would presumably be linked to the presence of 
a rotating magnetosphere in accreting neutron stars, which 
could directly imprint the spin frequency on the oscillations 
of the disk. 

The stellar spin frequency f * may couple to another pe- 
riodic motion already present in the system in more than one 
way, easily leading to the appearance of a beat frequency 
U2 = vi ± Vtf. However, if non-linearities are present, one- 
half of the spin frequency may appear in the data. We are 



fortunate in that nature has provided us with transient sys- 
tems where the accreting neutron star displays a coherent fre- 
quency, surely equal to its spin frequency. In the first of these 
systems to be closely studied, SAX J 1808, two kHz QPOs 
have been revealed, whose frequencies satisfy U2 — ui = 
z/*/2. In our view, this is a clear signature of non-Unearities 
in the system (Kluzniak et al. 2004). 

We may even try to give a specific model for the non- 
linear interaction. Elsewhere, we have shown that certain 
fluid configurations in the accretion flow, e.g. tori, can exe- 
cute two quasi-rigid oscillatory motions accurately described 
by equations nearly identical to that of a test particle in a cir- 
cular orbit, with the mode frequencies equal to the epicyclic 
frequencies in a certain orbit (Kluzniak and Abramowicz 
2002, Lee et al. 2004). Let us focus on one of these oscilla- 
tions, with the radial epicyclic frequency ujr- In the presence 
of the neutron star driver, it could be described by the equa- 
tion of a non-linear forced oscillator, with a forcing function 
f{t) that is periodic at the spin frequency, z^*: 

f + ujlr + a{z)r'^=f{t). (1) 

The configuration is also executing a vertical motion, 
whereby the value of the vertical co-ordinate z in this equa- 
tion oscillates at the frequency w^. We identify the observed 
frequencies with the frequencies of the oscillators: u\ = 
a;r/(27r) and V2 = ^^z/CStt). Reflection symmetry in the 
equatorial plane requires the expansion of the a{z) coeffi- 
cient to only have even powers in z. The leading coupUng 
term, then, has the form a2Z^r^, with a2a constant. Substitut- 
ing the harmonic time dependence of z, r, and / at frequen- 
cies respectively, f2, 1'l, and z/*, we see that the condition for 
forced resonance is 2(z/2 ± i^i) = i^*. For realistic models of 
the spinning neutron star in J 1808, the option 1^2 — 1^1 = i^*/2 
is viable. Some other possible forms of the non-Unear equa- 
tion also yield the same result. 

4. A model for the 3:2 frequency ratio 

The same oscillation model, even without the forcing term 
and the non-Unearity, allows a simple explanation for 
the 3:2 frequency ratio that is observed in the high fre- 
quency QPOs in black holes, as was noted in Abramowicz 
and Kluzniak 2001, Remillard et al. 2002, and as we argue, 
is also present in neutron stars. The epicyclic motions of the 
fluid configurations, e.g., of a torus are given by the equations 

f + ujlr = 0, (2) 

z + io,{rfz = 0. (3) 

In view of reflection symmetry in the z = plane, we 
have neglected the z variation of the radial epicyclic eigen- 
frequency. If only one oscillation is present, the frequencies 
are taken at a fixed equilibrium point, so that the solution 
is a simple harmonic oscillator motion. Now imagine that 
both oscillations are present. The radial motion induces har- 
monic variations of the eigenfrequency ujz, occuring at the 
frequency iVr, and the vertical motions can be resonantly ex- 
cited by the horizontal osciUations. This is a situation very 
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familiar to any child exciting a swing, formally the equation 
is known as the Mathieu equation, and the resonance (known 
as the parametric resonance) occurs when 

2 

UJr = -t^z, (4) 

n 

with rt = 1 , 2 , 3 ... an integer In general relativity, the first and 
strongest resonance occurs for n = 3, leading to a 3:2 ratio 
of the eigenfrequencies (Kluzniak and Abramowicz 2002), 
because, in additon to eq. (4), the relation 

must be satisfied. 

5. The Bursa plot 

The final question that must be discussed is whether neutron- 
star QPOs are indeed the same phenomenon as the black hole 
ones. This is not immediately obvious because the kHz fre- 
quencies in neutron stars vary. In fact the frequency ratio of 
these kHz QPOs departs from 3:2. There are three points to 
be made. 

One is that the histogram of frequency ratios for individ- 
ual neutron star sources, as well as for all of them jointly, 
peaks near the value 1.5 (Abramowicz, et al. 2003a). This is 
apparent in the inset of the figure (where actually the ratio of 
the lower frequency to the higher is plotted, so that the peak 
is close to 0.67). The black hole QPOs are much weaker so, 
at present, we only see the peak of the distribution. It is not 
know whether in black holes the weaker QPOs, presumably 
also present, have the same frequency, or whether their fre- 
quencies also depart from the 3:2 line. 

The second point is that the departure of the observed fre- 
quencies from the eigenfrequencies of the system, and of the 
ratio from the resonant 3:2 ratio, is a generic feature of non- 
linear resonance in coupled oscillators, so the line of cor- 
relation for an individual neutron star can be reproduced in 
numerical and analytic work (Abramowicz et al. 2003b, Re- 
busco 2004, Horak 2004). 

The third point, as is very apparent in the figure, is that 
the data points for at least some neutron stars seem to lie 
along straight lines on the frequency-frequency plot, with a 
common intersection point. This intersection point lies on the 
3:2 line and corresponds to frequencies which are higher than 
the black hole frequencies. This can be explained by assum- 
ing that these neutron stars have very similar masses, smaller 
than the (more varied) masses of the black holes. The lat- 
ter already show a l/Af frequency dependence discovered by 
Remillard and McClintock. The common intersection point 
of (some of) the neutron star lines greatly strengthens the ev- 
idence for the 3:2 resonance in neutron star QPOs, as well as 
for the strong gravitational field origin of the high frequency 
QPO phenomenon. 




100 200 300 400 500 600 700 800 900 1000 

V<iown [Hz] 

Fig. 1. The Bursa plot. The highest QPO frequency is plotted 
against its lower frequency counterpart. The black hole fre- 
quencies are from Remillard and McClintock 2005. The neu- 
tron star twin kHz QPO frequencies are based on published 
data (see Abramowicz et al. 2003a for references). Note that 
the neutron star points are scattered about a line intersecting 
the line of a 3:2 ratio in a single point ~ (600 Hz, 900 Hz), 
and that all black hole frequencies are lower, consistent with 
al/M dependence. The green line going through the Sco X- 
1 data has a slope of 3/4. This suggests a simple non-linear 
physical mechanism for its origin. 
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